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Abstract 



To establish the relationship between locomotory behavior and dy- 
namics of neural circuits in the nematode C. elegans we combined 
molecular and theoretical approaches. In particular, we quantitatively 
analyzed the motion of C. elegans with defective synaptic GABA and 
acetylcholine transmission, defective muscle calcium signaling, and de- 
fective muscles and cuticle structures, and compared the data with 
our systems level circuit model. The major experimental findings 
are: (i) anterior-to-posterior gradients of body bending flex for almost 
all strains both for forward and backward motion, and for neuronal 
mutants, also analogous weak gradients of undulatory frequency, (ii) 
existence of some form of neuromuscular (stretch receptor) feedback, 
(iii) invariance of neuromuscular wavelength, (iv) biphasic dependence 
of frequency on synaptic signaling, and (v) decrease of frequency with 
increase of the muscle time constant. Based on (i) we hypothesize that 
the Central Pattern Generator (CPG) is located in the head both for 
forward and backward motion. Points (i) and (ii) arc the starting 
assumptions for our theoretical model, whose dynamical patterns are 
qualitatively insensitive to the details of the CPG design if stretch re- 
ceptor feedback is sufficiently strong and slow. The model reveals that 
stretch receptor coupling in the body wall is critical for generation of 
the neuromuscular wave. Our model agrees with our behavioral data 
(iii), (iv), and (v), and with other pertinent published data, e.g., that 
frequency is an increasing function of muscle gap-junction coupling. 
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Introduction 



Caenorhabditis elegans nematode worms, with a small nervous system comprising 
only 302 neurons (White et al, 1986), move by generating an oscillatory neuromuscu- 
lar wave that alternates dorsal and ventral muscles (Brenner, 1974; Chalfie et al, 1985; 
Karbowski et al, 2006). The molecular, cellular, and network mechanisms of this oscilla- 
tory spatio-temporal activity are virtually unknown. Their understanding may provide 
insight about the relationship between neuromuscular dynamics and how behavior is 
created in these extensively genetically studied animals (Bargmann, 1998; Robert, 2003; 
de Bono and Maricq, 2005; Gray et al, 2005), and might be potentially relevant for other 
locomotory systems. It has proven difficult to address these questions using standard 
electrophysiological techniques in C. elegans because of its small neural sizes (Francis 
et al, 2003). In this paper, we use a combination of genetic perturbations, behavioral 
assays combined with a quantitative tracking system, and mathematical modeling to 
decipher dynamical properties of a detailed neuromuscular circuit relevant for C. elegans 
movement. By combining these approaches and building a mathematical circuit model 
we seek to bridge the gap between molecular/cellular and systems level understandings 
of undulatory locomotion (Fig. 1). 

We investigate specific questions related to the mechanisms that control body undu- 
lations and coordination: (i) How do different elements in the C. elegans neuromuscular 
circuit interact to produce both body oscillations and neuromuscular wave? In partic- 
ular, where is the primary oscillatory signal generated? (ii) How does mechanosen- 
sory feedback (stretch receptor coupling) affect locomotion? Does it play any role in 
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generating a neuromuscular wave? If so, how does its strength affect the wavelength 
(intermuscle phase lag) of muscle contractions? (iii) How does synaptic coupling be- 
tween neurons affect locomotory rhythm? Is there any qualitative difference between 
excitation and inhibition on oscillatory frequency? Are there optimal values for the 
synaptic couplings? (iv) How does gap-junction coupling between body-wall muscles 
and structural defects in muscles and cuticle affect movement? 

To address these questions we employed a parallel approach of collecting experi- 
mental data and making model predictions in an iterative manner. In particular, we 
quantitatively analyzed the motion of different neuronal and non-neuronal mutants by 
measuring their kinematic parameters (Fig. IB) and relating them to the dynamic 
properties of our circuit model. Investigated mutants included worms with decreased 
and increased synaptic GABA transmission, altered acetylcholine transmission, altered 
levels of calcium signaling in muscles, and worms with structural defects in muscles 
(myosin) and cuticle. 
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Materials and Methods 

Experimental pcirt 

Description of analyzed mutants. We examined the locomotion of several C. elegans 
mutants with affected inhibitory synaptic transmission (GABA), excitatory synaptic 
transmission (acetylcholine), and muscular function (calcium channels and myosin). 
We investigated mutants with both decreased and increased GABA function. For 
mutations effectively decreasing inhibition we studied: unc-25(el56), which encodes 
glutamic acid decarboxylase (GAD) - the biosynthetic enzyme for GABA production 
(Jin et al, 1999); unc-4.6(el77) {him-5(e 14.90) was in the background), which presum- 
ably plays a modulatory role in GABA packaging into vesicles (Schuske et al, 2004); 
and unc-18(e81) mutants carrying the syEx995[cho-l ::unc-18::yfpj transgene. unc-18 
functions as a facilitator of vesicles docking at presynaptic neurons (Weimer et al, 
2003) and is expressed in both cholinergic and GABAergic neurons (Gengyo-Ando et 
al, 1993). Therefore, in the unc-18(e81) strain carrying syEx995[cho-l::unc-18::yfp] 
we have restored unc-18 function specifically to cholinergic neurons using the cho-1 
promoter (Okuda et al, 2000), leaving the GABAergic neurons non-functional. For a 
mutation effectively increasing inhibition we studied slo-l(jsll8) mutants with extra- 
chromosomal Pacr-2::slo-l(-l-) transgene. SLO-1 encodes a Ca^"*" activated K"*" channel 
that presumably counteracts quantal synaptic neurotransmission both in excitatory and 
inhibitory neurons (Wang et al, 2001), and removal of this channel (in slo-1 mutants) 
leads to an increased neurotransmission above wild-type level. slo-l(jsll8) worms with 
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the Pacr-2::slo-l(+) transgenes (syEx996, syEx988, and syEx991) have elevated neu- 
rotransmitter release in inhibitory neurons, but they have restored wild-type release 
levels in excitatory neurons because Pacr-2::slo-l(+) drives expression of wild-type slo- 
1 in these neurons (Davies et al, 2003). We also studied three types of strains with 
presumably elevated acetylcholine (excitatory) synaptic signaling and wild-type level of 
GABA: slo-l(jsll8); Punc-25::slo-l, slo-l(jsll8); Punc-17::slo-l, and slo-l(jsll8); 
Pcho-l::slo-l. 

For mutations affecting muscular calcium function, we studied unc-68 mutants, as 
well as egl-19(n582) rescued with an unc-119::egl-19::yfp transgene. unc-68 gene en- 
codes ryanodine receptor channels that gate the release of calcium ions from internal 
stores in both body-wall muscle cells and in pharyngeal muscle (Maryon et al, 1996; 
Maryon et al, 1998). We analyzed the locomotion of loss-of-function unc-68(rll58) 
worms, which exhibit incomplete flaccid paralysis. The egl-19 gene plays an important 
role in regulating muscle excitation and encodes the al subunit of a homologue of verte- 
brate L-type voltage-activated Ca^"*" channels (Lee et al, 1997). Similarly, we analyzed 
the locomotion of egl(n582) mutants carrying an unc-119::egl-19::yfp transgene, which 
restores e^'Z-i 5 function specifically in neurons using the 'unc-ii5 promoter to drive pan- 
neuronal expression of egl-19 (Maduro and Pilgrim, 1995). This partial loss-of-function 
strain have slow muscle depolarization and feeble contractions (Lee et al, 1997; Jospin 
et al, 2002). 

For mutations affecting muscle and cuticle structure we analyzed unc-54, sqt-1, 
and BE109 worms, unc-54 worms carry mutations in the head region of body-wall 
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muscle myosin. We analyzed the following strains: RW130 unc-54(stl30), RW132 
unc-54(stl32), RW134 unc-54(stl34), RW135 unc-54(stl35), RW5008 unc-54(s95), and 
BC347 unc-54(s74 )■ These mutations are hypothesized to alter the contraction- relaxation 
cycle of the myosin-actin crossbridge formation by increasing its duration (Moerman 
and Fire, 1997). sqt-1 gene encodes cuticle collagen and we analyzed two strains: BElOl 
sqt-l(sclOl) and BE103 sqt-l(scl03). BE109 mutants have defective cuticle struts (J. 
Kramer, pers. communication). 

Construction of plasmids and strains is described in the Supplementary Information. 

Locomotory data extraction and processing. We video recorded and digitized the motion 
of young adult hermaphrodite worms (mutants and wild-type) 15-20 hr post mid-L4 de- 
velopmental stage on an agar plate with thin film of E. coli OP50 bacteria and LB media. 
The video recording and data extraction was done using a device specially designed for 
studying Caenorhabditis locomotion (Cronin et al, 2005). We collected 5 minutes of 
video per worm, extracting digital locomotion data from the middle 4 minutes. Such 
4 minute windows average over many possible sensory influences that can vary among 
worms and thus statistically minimize the variability of external conditions. From these 
data we derived distributions of the undulatory frequency and bending flex by dividing 
worm's body into 12 sections which served as worm's system of coordinates. Local flex 
angle was determined as described in Fig. IB. To derive the bending frequencies for 
each articulation point we flrst calculate a time-ordered matrix of angles between each 
pair of adjacent body sections (that is, at each articulation point). Next we apply Mat- 
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lab's spectrogram function (specgram) to the vector of changing bend angles at each 
articulation point, which yields a time-ordered matrix of the relative energies of the 
changing bend angles' component frequencies. The highest magnitude, non-constant 
component frequency from each time window is that window's characteristic bending 
frequency. Frequencies and flexes for each articulation point in Figs. 3-6 are population 
averages. In Figs. 5 and 6 we used non-parametric sign test for paired samples to 
determine if average frequencies between different mutants are significantly different. 

Theoretical peirt 

Circuit model. All known locomotory systems use some form of a Central Pattern Gen- 
erator (CPG) circuit to generate rhythm (Grillner, 1975; Delcomyn, 1980; Marder and 
Calabrese, 1996; Nusbaum and Beenhakker, 2002). Mechanosensory feedback usually 
alters the rhythm; however, its importance varies in different animals and even in the 
same animal under different conditions (Marder et al, 2005). At one extreme, it can 
play a minor modulatory role (Grillner, 1975; Delcomyn, 1980), at the other, it can be 
critical for body coordination (Friesen and Gang, 2001; Akay et al, 2004). 

The global circuit responsible for C. elegans forward locomotion is constructed based 
on neuroanatomical data (Ghalfie et al, 1985; White et al, 1986). It is composed of the 
head and body wall neural networks (Fig. 2A,B) connected partly via two interneurons 
AVE and PVG, which are implicated as controllers of forward movement (Ghalfie et al, 
1985), and partly via muscle gap junctions. AVE and PVG receive massive synaptic 
input from the head neurons and target only excitatory motor neurons (type E neurons: 
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ventral VB and dorsal DB) in the body wall (White et al, 1986). AVB connects them 
via gap junctions and PVC via chemical synapses. We assume, following theoretical 
predictions (Wicks et al, 1996), that the latter synapses are excitatory. 

The real head neural network is composed of about 200 neurons with mostly un- 
known polarities and complicated wiring diagram (White el al, 1986). Inclusion of all 
of these elements in our circuit model would increase its complexity (the number of 
parameters) by at least an order of magnitude, which would obscure the basic mech- 
anisms. Instead, we decided to depict the head network in a simplified manner with 
lesser but essential number of elements. This simplified approach has an advantage of 
being conceptually comprehensible. Thus our head neural network is comprised of head 
interneurons and ventral and dorsal head motor neurons (Fig. 2B). The head motor 
neurons contain both excitatory neurons (White et al, 1986) and a pair of GABA-ergic 
inhibitory RME neurons (Schuske et al, 2004). Because connectivity patterns between 
actual head interneurons form numerous feedback loops (e.g., AIZ AIA ^ AWA — >■ 
AIZ and RIB ^ RIG ^ URY RIB; see White et al, 1986), we mimicked this phe- 
nomenon by two (dorsal and ventral) interneuron loops X— > Y— > Z— > X. These loops 
can in principle generate a primary rhythm in the head, which is then imposed on 
downstream networks along the body via muscle coupling and/or AVB and PVC in- 
terneurons. The loop X, Y, Z is the minimal circuit in the C. elegans head interneurons 
that can generate robust oscillations. We include mechanosensory feedback in our head 
network model following suggestions that a pair of head interneurons SAA may contain 
stretch receptors on their long dendrites (White et al, 1986). 
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The body wall forward motion network (White et al, 1986) is composed of ventral 
and dorsal sub-circuits that contain motor neurons, body wall longitudinal muscles, and 
feedback loops associated with mechanosensation (Fig. 2A). The neuronal part of each 
sub-circuit consists of several cholinergic excitatory B motor neurons (ventral VB and 
dorsal DB), and several GABA-ergic inhibitory motor neurons (type D neurons: ventral 
VD and dorsal DD) (Chalfie et al, 1985; Mclntire et al, 1993; Schuske et al, 2004). 
Apart from gap junctions from AVB and chemical synapses from PVC, the excitatory 
B neurons receive synaptic input from DVA interneuron (DVA synapses on some B 
neurons, see White et al, 1986). Anatomically, the role of this interneuron in locomotion 
is not clear, and it was suggested (Li et al, 2006) that it might serve as a neuromodulator 
involved in detecting body stretch, with no definite synaptic polarity. We assume that 
the overall influence of this neuron is inhibitory. This inhibitory input is necessary to 
prevent spontaneous oscillations in the body-wall circuit when it is disconnected from 
the head. Coupling between ventral and dorsal motor sub-circuits is mediated partly via 
excitatory synapses terminating on inhibitory D motor neurons of the complementary 
circuit (Fig. 2A), and partly via AVB gap-junctions. An effective cross-inhibition 
(B D) between the two sub-circuits amplifies their alternating oscillatory activities 
(Fig. 2A, and below). Neighboring motor neurons, as well as neighboring muscles are 
connected via gap junctions (D. Hall, pers. commun). Mechanosensory feedback in the 
body wall, similarly as in the head, is associated with hypothetical stretch receptors. 
Following the Byerly and Russell hypothesis (in Chalfie and White, 1988), we assume 
that stretch receptors are located on the extended dendritic processes of the excitatory 
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B motor neurons (VB and DB). In the forward motion circuit, all these dendrites are 
directed posteriorly. We assume that mechanosensors transduce local body stretch 
along excitatory motor neuron dendrites to the somas of these neurons by polarizing 
them, since mechano-stimulation is known to increase cell excitability (Goodman and 
Schwarz, 2003). Thus, mechanosensory feedback coupling along the body wall is non- 
local, excitatory, and unidirectional in nature. 

The theoretical part of this study deals almost exclusively with forward motion, 
which is a dominant behavior in C. elegans. Backward motion is considered only briefly 
in the context of neuromuscular wave directionality (see below) . The backward motion 
is controlled by a complementary system of 3 interneurons (AVA, AVD, AVE) and type 
A excitatory motor neurons (ventral VA and dorsal DA). The latter neurons have long 
dendrites directed anteriorly, i.e., opposite to the corresponding motor neurons in the 
forward motion circuit (Chalfie et al, 1985). The behavioral transitions between forward 
and backward motion will be addressed in a future study. 

Equations used in the theoretical model. We assume that the activities of motor neu- 
rons are graded, since experimental data show the absence of action potentials in the 
sensory neuron ASER (Goodman et al, 1998), which is consistent with the lack of 
voltage-activated sodium channels in the C. elegans genome (Bargmann, 1998). Ad- 
ditional support for this assumption comes from neurophysiological recordings in a 
related, larger, nematode Ascaris, which has a similar nervous system, and lacks ac- 
tion potentials in all recorded neurons (Davis and Stretton, 1989). Below we provide 
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equations governing activities of all neurons, muscles, and stretch receptors included in 
the model. These equations were solved on a Linux-system computer using a standard 
second-order Runge-Kutta method in Fortran (F77). 

The dynamics of head interneurons controlling the ventral side of the head are 
represented by: 



dX 



dYy 



dZ„ 



"^z-^ = —Zv + WzyHzyiYy) + Wzx^zx{Xd) + gavb{Vavb " Zy) , 



(2) 
(3) 



where Xy and Yy denote the activities of two excitatory ventral head interneurons, X^ 
denotes the activity of a corresponding head dorsal interneuron, Zy is the activity of 
ventral inhibitory interneuron, and Tx, Ty, Tz are corresponding time constants. Shv is 
the head stretch receptor activity situated on the ventral interneuron X^, and is 
the membrane potential of head inhibitory ventral motor neuron. Vavb is the membrane 
potential of AVB interneuron and gavb is the gap junction strength between AVB and X, 
Y, and Z. The parameter Ci denotes some input coming from other head neurons and 
it can be either excitatory (positive) or inhibitory (negative), depending on the worm's 
"internal state". The parameters characterize the strength of chemical synapses 
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between circuit elements from /3 to a, and most of the genetic mutations we examined 
affect precisely these parameters. The functions Hap are positively defined sigmoidal 
functions of the form: Hq,^(x) = (1 + tanh[(2; — Oai3)/r]ai3]), where is the threshold 
for activation and rja/s characterizes the steepness of the nonlinearity. The equations for 
the dorsal side of the head are analogous with a subscript substitution v ^ d. 

The dynamics of the ventral side head motor neurons, muscles, and stretch receptor 
feedback are given by: 



Te^^ = -Ehv + Weyii^yiYy), (4) 



^i-^ = -hv + WieUi^iEhv), (5) 
mm 

dt 

-WmiH.ini{Ihd) + QmiMy^l " M^^): (6) 



Ts^^ = -Shv + Wsm [Hsm(Mfed - Mhv) - 1] , (7) 



where E^v and Im are the membrane potentials of head ventral excitatory and dorsal 
inhibitory motor neurons respectively, Mhv and Mhd are head ventral and dorsal mus- 
cle membrane potentials, and Tg, Tj, Tm, Ts arc their respective time constants. The 
steady-state value of the stretch receptor activity Shv should be proportional to the 
degree of head bending, which we assume is a sigmoidal function of the differences in 
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ventral and dorsal muscle activities. Note that stretch receptors are activated if there 
is initial asymmetry in activities of ventral and dorsal muscles. Apart from synaptic 
input, a muscle cell is also activated by a voltage-activated calcium channels, which in 
the model is represented by a self-promoting HmmiMhv) term with the coupling Wmm 
characterizing the calcium signaling strength. Identical equations describing the dorsal 
side are obtained with a subscript substitution hv hd. 

The dynamics of the two forward motion interneurons AVB and PVC are given by: 



'^avb ^avb + '^avb,x [Havb,x 

N 

~l~'^aii6,piicHavb,pvc(^i;c) ^" 9avb,e ^ ' ~l~ 2V^j,j,) 

i=l 

+gavb (X^ + Y^ + Z^+Xd + Yd + Zd- GVavb) , (8) 



TpVC ^pVC ~l~ '^pVC,X [Hpvc,x(^l) ~l~ ~l~ Hpvc,x(A^d ~l~ ^^)] ) (9) 



where Vavb and Vp^c are the membrane potentials of AVB and PVC interneurons and 
Tavb,Tpvc are their respective time constants. Ey^i and E^^i denote the membrane po- 
tentials of the ith ventral (VB type) and dorsal (DB type) excitatory motor neurons, 
respectively. The parameter gavb,e is the gap-junction coupling between AVB and the 
body wall excitatory motor neurons. N denotes the numbers of excitatory and in- 
hibitory motor neurons on both sides of the body wall; in simulations we take N ~ 8. 
The ventral side of the body wall nervous system is represented by a set of equations: 
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) , (10) 



~ ~-^v,i + Wie^ie{Ed,i) + Qi {Iv,i+1 + Iv,i-1 — '2Iv,i) , (H) 

-«^miH„i(4,i) + Qm {M^,i+i + - 2M^^i) , (12) 



dt 



-Sy,i + Wsm [iism{Md,i - M^^i) - 1] , (13) 



where I^^i is the membrane potential of the ith ventral inhibitory motor neuron (VD 
type), and Sy^i is the stretch receptor activity related, as in the head, to the relative 
activities of the ventral M^, j and dorsal M^^i muscles. By the nature of connectivity in 
our network, My^o = ^hv and M^^ = M^- Note that the ith excitatory motor neuron 
Ey^i (corresponding to VB type) is affected by the stretch receptors at (i + l)th position 
in the forward motion circuit. (In the backward motion circuit E^^i corresponds to VA 
neuron and is affected by the stretch receptor at {i — l)th position.) The parameters g^i 
Qi, and Qm are the gap-junction couplings among excitatory motor neurons, inhibitory 
motor neurons, and among muscles, respectively. The parameter C2 denotes inhibition 
coming to the excitatory motor neurons (VB, DB) from the DVA interneuron, and 
its negative value prevents spontaneuos oscillations in the body wall circuit when it 
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is disconnected from the head circuit. These spontaneous oscillations emerge because 
of the numerous closed loops in the body wall, which contains both excitatory and 
inhibitory elements (VB DD ^ Ma ^ Sa ^ DB ^ YD ^ ^ ^ VB), see 
Fig. 2A. Identical equations describe a dorsal part of the body wall circuit, with a 
substitution v d. The ventral inhibitory motor neuron j (VD type) is activated by 
the activity of the dorsal excitatory motor neuron E^^i (DB type). 

In simulations we rescale uniformly almost all synapses both in the head and in the 
body wall by two factors: q^.^ corresponding to excitatory synapses, and qin correspond- 
ing to inhibitory GABA synapses. Both parameters are variable and characterize global 
levels of excitation and inhibition, respectively. The values of parameters used in simu- 
lations arc presented in Table 1. Values of these parameters were chosen arbitrarily and 
a global circuit function (i.e. sustained oscillations), to a large extent, does not depend 
on their precise values. Oscillations in the circuit emerge for c\ between 0.5 and 2.5. 
Models A and B differ only in one aspect, i.e., the degree of nonlinearity in interactions 
(in functions H) between head interneurons X, Z. Model A was obtained by tak- 
ing the corresponding steepness parameters r] to be small (high nonlinearity), whereas 
model B was obtained by increasing values of these parameters (low nonlinearity). The 
rest of the parameters are the same for both models. In all theoretical figures we used 
the parameters from Table 1, unless indicated otherwise. 

The Fortran code for simulating of Eqs. (1)-(13) will be deposited in ModelDB 
database. 
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Results 



To begin an analysis of how the circuit controls locomotion, we chose specific genetic 
perturbations that affect easily identifiable classes of parameters in the neuroanatomical 
circuit. These parameters can be manipulated in our theoretical model, enabling a 
comparison between theory and experiment. 

Experimental results 

The major experimental findings detailed below are: (i) anterior-to-posterior spatial 
gradients of flex angle for almost all strains, and for neuronal mutants, also analo- 
gous weak gradients of frequency, both for forward and backward motion, (ii) cuticle 
and muscle (non-neuronal) mutants affect locomotion characteristics, suggesting the 
existence of some form of mechanosensory feedback, (iii) invariance of wavelength char- 
acterizing the neuromuscular wave, (iv) biphasic dependence of undulatory frequency 
on synaptic GABA signaling, and a plausible biphasic dependence on acetylcholine sig- 
naling, and (v) strong decrease of undulatory frequency with increasing muscle time 
constant characterizing the cycle of myosin-actin crossbridge formation. 

Spatial distribution of undulatory frequency and bending fiex. The spatial distribution 
of some locomotory parameters may provide clues about how the motion is generated 
and coordinated in C. elegans worms. We measured frequency of undulations and 
bending flex, deflned as a maximal local bending angle of worm's body, as a function 
of their distance from the head. The frequency for wild-type and most non-neuronal 
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mutants is either constant or nearly constant both for forward and backward motion 
(Fig. 3A,B, and Supplementary Table Tl), with the exception of sqt-1 (lacking one 
cuticle collagen), for which the frequency decreases significantly toward the tail (e.g., in 
sqt-1 (scl03), the head and tail oscillate at respectively 0.29 ± 0.07 Hz and 0.22 ± 0.09 
Hz in forward motion, and corresponding values in backward motion are: 0.22 ± 0.05 
Hz and 0.12 ±0.04 Hz; see Suppl. Table Tl). In contrast, all neuronal mutants display 
a monotonic decay of the frequency directed posteriorly, in some cases up to 50 % (see 
Suppl. Table Tl). Specifically, we obtained the following frequencies in the head and 
tail for neuronal mutants. For unc-25(el56): 0.36 ±0.07 Hz in the head and 0.23 ±0.06 
Hz in the tail (forward motion), and 0.20±0.08 Hz in the head and 0.13±0.03 Hz in the 
tail (backward motion); for unc-46(el77): 0.29 ± 0.03 Hz in the head and 0.16 ± 0.07 
Hz in the tail (forward motion), and 0.26 ± 0.02 Hz in the head and 0.13 ± 0.04 Hz in 
the tail (backward motion); for unc-18(e81); cho-l::unc-18:yfp: 0.20 ± 0.08 Hz in the 
head and 0.09 ± 0.04 Hz in the tail (forward motion), and 0.16 ± 0.05 Hz in the head 
and 0.06 ±0.02 Hz in the tail (backward motion); for slo-l(jsll8): 0.38 ±0.07 Hz in the 
head and 0.29 ± 0.08 Hz in the tail (forward motion), and 0.32 ± 0.07 Hz in the head 
and 0.16 ±0.06 Hz in the tail (backward motion). These data indicate that for neuronal 
mutants tail oscillations have substantially smaller frequencies than head oscillations. 

The behavior of the fiex is more universal (Fig. 3C,D and Suppl. Table Tl). For 
almost all mutants and wild-type worms, and for both forward and backward motion, 
there is a pronounced decay of the fiex from the head up to some point close to the 
tail. From that point the fiex either stays constant or shghtly increases towards the tail. 
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Overall, the flex and position exhibit statistically signiflcant strong negative correlations 
(Fig. 3 C,D). Thus, surprisingly, the opposite motions yield similarly directed anterior- 
to-posterior gradients of the flex. 

Some non-neuronal mutants dramatically alter frequency. We flnd that it is relatively 
easy to change the undulatory frequency of worms by genetically modifying some of non- 
neuronal parameters. As an example, mutants with affected cuticle function (collagen 
defective sqt-1 mutants or strut defective BE109 mutants) and mutants with altered 
muscle function (myosin defective unc-54 mutants) exhibit different locomotory outputs 
than the wild- type worms, with substantially reduced frequency (Fig. 4). Also, results 
(see below) on mutants with defects in muscle calcium signaling {unc-68) show that their 
frequency is significantly lower than wild-type. The fact that non-neuronal mutants have 
profoundly altered locomotory parameters, in particular frequency and flex, indicates 
that there must be some feedback coupling between the nervous system of C. elegans 
and its muscular and cutico-skeletal systems. 

Conservation of wavelength. We measured the wavelength of neuromuscular wave in 
wild type and several mutants (Table 2). The ratio of the wavelength to the actual 
worm's body length is relatively constant across different worms and appears to be 
identical for forward and backward motion. Since these worms move with different 
frequencies, the data in Table 2 suggest that wavelength is essentially frequency inde- 
pendent. 
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Biphasic dependence of undulatory frequency on GAB A (inhibitory) transmission. We 
analyzed the motion of mutants defective in GABA function primarily in the ven- 
tral cord motor neurons. All three uncoordinated neuronal mutants - unc-25(el56), 
unc-46(el77), and unc-18(e81) rescued with transgene syEx98 l[cho-l::unc- 18: :yfp] in 
cholinergic neurons - have decreased GABA transmission as compared to wild type 
(Mclntire et al, 1993; Weimer et al, 2003). These animals have also much lower fre- 
quencies of undulations than wild-type worms (Fig. 5 A and Suppl. Table Tl, p < 0.05). 

We also studied the behavior of worms with increased GABA transmission. Since 
no known mutants have this property exclusively, we constructed and analyzed sev- 
eral transgenic strains with altered GABAergic and cholinergic signaling based on slo- 
l(jsll8) mutation, which elevates globally both types of neurotransmission by eliminat- 
ing a particular calcium activated potassium channel (Wang et al, 2001) (Suppl. Table 
T2). By using different promoters we constructed worms that restore either excitation 
or inhibition to normal levels by making channel functional in either cholinergic 
or GABA-ergic neurons. The slo-l(jsll8); Pacr-2::slo-l worms have presumably res- 
cued cholinergic transmission to the wild-type level, but GABA transmission is still 
elevated. These worms are significantly slower {p < 0.05) than wild-type worms (Fig. 
5B). For consistency, we also compared slo-l(jsll8); Pacr-2::slo-l animals to two strains 
with presumably wild-type levels of both GABA and acetylchohne, i.e., slo-l(jsll8); 
Pacr-2::slo-l; Punc-25::slo-l and slo-l(jsll8); Psnb-l::slo-l (Fig. 5C). In this case 
slo-l(jsll8); Pacr-2::slo-l are also significantly slower [p < 0.05). Thus, the behavior 
of worms with increased GABA signaling (Figs. 5B,C) is qualitatively similar to the 
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behavior of worms with decreased GABA (Fig. 5A): the undulatory frequency of worms 
with altered GABA transmission, either below or above the wild-type level, is always 
smaller than wild type. 

To summarize, all these results clearly suggest a non-monotonic (biphasic) depen- 
dence of the undulatory frequency on the inhibitory synaptic GABA signaling in the 
worm. 

Dependence of undulatory frequency on acetylcholine (excitatory) transmission. Various 
constructs with slo-1 enable us to determine the dependence of frequency on excitatory 
coupling in C. elegans (Suppl. Table T2). slo-l(jsll8); Pacr-2::slo-l animals with pre- 
sumably rescued levels of acetylcholine transmission to the wild- type level but elevated 
GABA transmission are either slightly faster or slightly slower than the control slo- 
l(jsll8); Pacr-2::gfp with both increased GABA and increased acetylcholine signaling 
(Fig. 6A). Thus, a decrease of excitatory transmission does not lead to a clear-cut phe- 
notype, which might suggest a possible non-monotonic dependence of the frequency on 
a global synaptic excitation. In the case of worms with increased acetylcholine signaling 
and presumably wild-type level of GABA (three strains: slo-1 (jsl 18); Punc-25::slo-l, 
and slo-l(jsll8); Punc-17::slo-l, and slo-l(jsll8); Pcho-l::slo-l), we observe their 
slightly lower frequencies than wild type (Fig. 6B). This case, in turn, suggests that 
increasing global excitation leads to a slight decrease in undulatory frequency (see, 
however, "Comparison of the theory with the experimental data" section below). 

Frequency decreases sharply with increasing muscle time constant. Fig. 3 and Supple- 
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mental Table Tl show that all strains of unc-54 mutants move with frequencies much 
lower than wild type. These mutants have affected muscle myosin in the body wall, 
which presumably increases the time constant controlling contraction- relaxation cycle of 
myosin-actin crossbridge. Thus, increase in muscle time constant leads to a pronounced 
decrease in the undulatory frequency both for forward and backward motion. 

Dependence of undulatory frequency on muscle calcium signaling. We studied locomo- 
tion of C. elegans mutants with defective calcium signaling in body wall muscle cells. 
Null unc-68 mutants, which have defective ryanodine receptor channels responsible for 
release of stored calcium ions, display abnormal locomotion (Maryon et al, 1996; Maryon 
et al, 1998). We analyzed the motion of loss-of-function unc-68(rll58) worms with re- 
duced levels of calcium in muscles, and found that their frequencies of undulations are 
smaller than that of the wild- type worms (Suppl. Table Tl and Fig. SI), in agreement 
with Maryon et al. (1996, 1998). 

We also analyzed locomotion of loss-of-function egl-19(n582) animals carrying a 
transgene expressing EGL-19-YFP in neurons {unc-119::egl-19::yfp). This strain has a 
defective voltage-gated Ca+^ channel in body wall muscles but has its function restored 
in neurons. These worms have less muscle calcium influx (Jospin et al, 2002) and lower 
frequency of undulations than wild type (Suppl. Table Tl and Fig. SI), similarly to 
the unc-68 worms. 

These data suggest that decreasing calcium signahng in muscle cells leads to decrease 
in the frequency of body oscillations. Moreover, because unc-68 is expressed almost 
exclusively in muscles (Maryon et al, 1998), and egl-19(n582); unc-119::egl-19::yfp has a 
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defect also only in muscles, these results provide an additional support for our hypothesis 
about the existence of a feedback coupling between body wall muscles and the neural 
circuits. 

Theoretical results 

The major theoretical findings detailed below are: (i) the dynamic properties of the 
global circuit practically do not depend on CPG design in the head if stretch receptor 
feedback is sufficiently strong and slow, (ii) the neuromuscular wave is generated due 
to a nonlocal stretch receptor feedback coupling in the body wall, and its wavelength is 
almost frequency and stretch receptor coupling (if strong enough) independent, (iii) the 
frequency of oscillations depends in a biphasic manner on both inhibitory and excitatory 
synaptic coupling, (iv) the frequency increases monotonically with muscle gap-junction 
coupling strength, and (v) the frequency decreases with increasing time constants of the 
circuit elements. We also predict that frequency should depend non-monotonically on 
muscle calcium signaling. We made no attempt to fit quantitatively the experimental 
data, since many physiological parameters in C. elegans are unknown. In this sense, 
our approach is qualitative. However, we adjusted model parameters, especially r^, to 
obtain frequencies of the same order of magnitude as in experiments. 

Oscillation generation. We assume that primary oscillations necessary for a regular 
movement (i.e., spontaneous and not caused by an outside mechanical stimulation) 
are generated in the head. The arguments, based on the anterior-to-posterior gradient 
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of the flex (Fig. 3), supporting the head as a hkely primary CPG are presented in 
the Discussion. We consider two extreme design versions of the head CPG. In model 
A, the stretch receptor feedback plays a minor modulatory role in affecting neurally 
generated rhythm (neural component dominates). On the other hand, in model B, the 
stretch receptor feedback plays a critical role in rhythm generation (reflex component 
dominates). In model A oscillations are generated in the head interneuron circuit, i.e., 
in the loop X — >^ Y — >■ Z — > X, provided nonlinearities in the interactions between 
these interneurons are sufficiently strong (corresponding sigmoidal functions Ha/d have 
sharp thresholds for activation). The emerging oscillatory patterns are presented in 
Fig. 7A,C. One of the roles of the stretch receptor feedback in model A is to slow 
down the emerging rhythm as compared to a purely neural oscillator by introducing 
temporal delays. (This is because the characteristic membrane time constant of C. 
elegans neurons should be of the order of 100 msec, as it is in morphologically similar 
Ascaris neurons (Davis and Stretton, 1989) and much larger mammalian neurons. As 
a result, a purely neural rhythm would give ~ 10 Hz body undulations, which is 20-30 
times larger than that observed experimentally; Suppl. Table Tl.) In model B, the 
primary oscillations in the head arc generated because of the presence of the stretch 
receptor feedback and its interaction with neurons, i.e. in the loop X^Y— >E— >M— > 
S — >^ X. Oscillations in this model are generated even when the interneuron circuit (X, Y, 
Z) is intrinsically non-oscillatory, i.e., when interactions between interneurons are only 
weakly nonlinear or linear. In this case, however, the prerequisite for the appearance 
of oscillations is initial asymmetry in activities between ventral and dorsal sides, which 
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is necessary for engaging the stretch receptor feedback. Such activity asymmetry can 
be faciUtated by nonequal number of dorsal and ventral motor neurons present in C. 
elegans. The emerging oscillatory patterns in model B are presented in Fig. 7B,D. 

In both models of CPG, the crucial element for oscillation generation is a net non- 
linear negative feedback loop destabilizing the steady state and stabilizing a limit cycle 
(periodic state). The frequency of the emerging oscillations depends on neuronal and 
non- neuronal parameters characterizing the oscillatory circuit. The constant input ci 
(see Eq. (1) in Materials and Methods) coming to the head interneurons from other 
head neurons or/and sensory input can abolish oscillations if it is either too inhibitory 
or too excitatory. Thus, depending on the activity level of the head neural networks, 
a worm can decide when to start or to stop moving. Our theoretical analysis (below) 
shows that models A and B give very similar results for all model parameters, except 
for the coupling strength and time constant of the stretch receptors. 

Oscillations generated in the head are imposed on the body wall nervous system 
primarily through interconnections of the body wall muscles with the head muscles, 
and then via body wall stretch receptors and synapses. However, the electric coupling 
between muscles and synaptic coupling between neurons cannot be too weak, since 
then oscillations along the body wall are unsustainable and their amplitude (waveform) 
decays with a distance from the head (Fig. 7A,B, lower panels). In our model the 
stability for the oscillations is provided by the interneurons AVE and PVC Since these 
neurons receive a massive oscillatory input with different phases in the head (from X, 
Y, Z) , their resulting activities are only weakly- or non-oscillatory, because of a partial 
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or complete cancellation of the input oscillatory phases. Nearly constant input from 
AVB and PVC to excitatory motor neurons overcomes their inhibition from the DVA 
neuron and stabilizes undulatory oscillations, with ventral and dorsal sides in anti-phase 
(see also Discussion). This input in combination with sufficiently strong synapses and 
gap-junctions ensure that the oscillatory waveform does not decay with a distance from 
the head. Consequently, the distribution of frequencies along the body is homogeneous 
despite the fact that rhythm is generated in the head. This homogeneity is consistent 
with the experimental results for wild- type and all non-neuronal mutants (Figs. 3, 4 
and 7; Suppl. Table Tl). For the distribution of frequency of neuronal mutants with 
weak synapses see the Discussion. 

Wave generation and wavelength constancy. Undulatory movement also requires a neu- 
romuscular wave (Karbowski et al, 2006). Our model for forward locomotion can gener- 
ate such a wave due to a non-local stretch receptor activity transduced to the body wall 
somas of excitatory motor neurons via their long unidirectional dendrites (Fig. 2). With 
a posterior stretch receptor coupling the wave propagates posteriorly (corresponding to 
forward motion), while with an anterior stretch receptor coupling the wave propagates 
anteriorly (corresponding to backward motion) (Fig. 8A,B). Thus, the directionality 
of the wave (and thus motion direction) depends only on the directionality of the long 
dendrites. Unidirectional stretch receptor feedback along the body wall is critical for lo- 
comotion, since when it is absent the wave does not propagate. In this case oscillations 
are also abolished in model B, although in model A they are still present (due to an 
intrinsically oscillatory interneuron loop X — > Y — > Z — > X). The correlation between 
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directions of the wave and a hypothetical stretch receptor couphng is consistent with 
the neuroanatomical details of the corresponding body wall excitatory motor neurons 
of type A and B, controlling backward and forward motion, respectively. These neurons 
have long undifferentiated dendrites directed in opposite directions (White et al, 1986). 
The phase differences between neighboring muscle activities along the body wall are to 
a large extent stable over time and essentially frequency independent (Fig. 8C). Since, 
in general, a phase lag is inversely proportional to wavelength, this implies that the 
wavelength is also frequency independent. The intermuscle phase lag depends, how- 
ever, on the stretch receptor coupling, although very weakly if coupling is sufficiently 
strong (Fig. 8D). 

Biphasic dependence of frequency on synaptic coupling. Chemical coupling via synapses 
is one of the critical parameters in the circuit, since it can potentially affect numer- 
ous kinematic characteristics. Our model makes predictions about the influence of 
inhibitory and excitatory transmission on frequency. In Fig. 9A,B, we plot the de- 
pendence of the frequency on the inhibitory synaptic strength at the neuromuscular 
junction. This dependence exhibits a non-monotonic (biphasic) relationship. For weak 
inhibitory strength, the frequency is small and it increases with increasing inhibition 
up to some maximal value, and then it decreases. 

With respect to excitation, our model shows that frequency is a biphasic function 
of the global excitatory synaptic strength (Fig. 9C,D), similar to the above case of 
inhibition. When the excitation level is either above or below some critical values, the 
oscillations become unstable and movement disappears. For the consistency of these 
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results with the experimental data see the section "Comparison of the theory with the 
experimental data". 

The origin of the biphasic relationships between the frequency and inhibitory and 
excitatory synaptic strengths lies in the shapes of temporal muscle voltage character- 
istics (Fig. 9E,F). For weak inhibition (or strong excitation) muscle voltage increases 
sharply but its relaxation phase is prolonged as compared to the case of intermediate 
levels of inhibition (or excitation). In contrast, for strong inhibition (or weak excitation) 
the muscle voltage increases slowly but its decay phase is much faster. In both regimes, 
corresponding longer decay or rise phases of the muscle voltage cause the frequency 
drop in comparison to the intermediate regime (Fig. 9E,F). 

Frequency vs. gap-junction coupling. Electrical coupling via gap junctions plays also 
a significant role in coordination of the circuit components. Models A and B both 
predict that frequency is a monotonically increasing function of the muscle gap junction 
conductance strength (Fig. 10A,B). Moreover, when this gap junction coupling is too 
weak, the oscillatory signal coming from the head fails to propagate down the body, 
i.e., the neuromuscular wave is attenuated, causing defective locomotion. On the other 
hand, when this coupling is too strong the wave disappears because the whole system 
is synchronized. Thus, for a robust wave the gap-junction coupling should be in an 
intermediate range. 

Frequency vs. muscle calcium signaling. In our model, the amount of calcium influx in 
muscles is controlled by the coupling Wmm (Eqs. 6 and 12). In Fig. 10C,D, we display 
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the biphasic or triphasic dependence of frequency on this amphtude. At some very 
high level of calcium signaling oscillations become unstable and the system converges 
to a steady-state, in which undulatory motion is abohshed. The curves in Fig. 10C,D 
are qualitatively distinct from the biphasic curves for global excitation and inhibition 
presented in Fig. 9, since the latter exhibit pronounced local maxima. 

Frequency vs. stretch receptor coupling. One of the roles of the stretch receptor feedback 
couphng is to modulate the undulatory rhythm. We find that models A and B differ 
only for a weak stretch receptor coupling. However, if the coupling is sufficiently strong, 
both models show that frequency increases slightly with increasing the coupling (Fig. 
11). 

Frequency vs. time constants. We also tested how the frequency of oscillations depends 
on various time constants associated with different dynamic variables. In general, an 
increase in any of the time constants causes decrease in frequency, with a variable 
sensitivity among different time constants. As an example, in Fig. 12 we show the 
dependence of frequency on the muscle membrane time constant (Fig. 12A), and on 
the stretch receptor time constant Tg (Fig. 12B). Models A and B differ only for small 
values of time constants. 

Compcirison of the theory with the experimental data 

Below we make a qualitative comparison of several of our theoretical results to either 
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our data or to the data already published. 

• Invariance of neuromuscular wavelength. This experimental result (Ta- 
ble 2; see Karbowski et al. (2006) for a larger data set) is consistent 
with our theory, which shows that intermuscle phase lag is practically 
frequency independent (Fig. 8C), and is only weakly stretch receptor 
coupling dependent if this coupling is sufficiently strong (Fig. 8D). The 
latter result can explain wavelength constancy across different nema- 
tode species (Karbowski et al, 2006). 

• Biphasic dependence of frequency on synaptic coupling. For inhibitory 
coupling this experimental result (Fig. 5) is consistent with our theory 
(Fig. 9A,B). Moreover, the theoretical curve suggests that wild-type 
worms may have GABA coupling parameters close to the optimal val- 
ues. 

In the case of excitatory coupling, data in Fig. 6A for acetylcholine 
signaling are not in conflict with the theoretical biphasic relationship 
(Fig. 9C,D), although they are not conclusive. However, the support 
for the biphasic relation between frequency and acetylcholine signaling 
comes from combining the experimental results in Fig. 6B, containing 
slo-1 worms with transgenes, with our previously published data (Kar- 
bowski et al, 2006) on hyperactive goa-l(nll34), goa-l(syl92), and 
egl-30(tg26) mutants (Segalat et al, 1995; Mendel et al, 1995). Both 
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of these groups have increased acetylchohne synaptic signahng above 
wild-type level but their corresponding undulatory frequencies show 
the opposite behavior. In the case of slo-1 worms with transgenes, the 
frequency slightly decreases as compared to animals with the wild-type 
levels of neurotransmission (Fig. 6B), while in the case of hyperactive 
mutants the frequency increases (Karbowski et al, 2006). This behav- 
ior clearly suggests a non-monotonic, possibly biphasic dependence, in 
agreement with our theoretical results (Fig. 9C,D). 

• Frequency increases (decreases) monotonically with increasing (decreas- 
ing) muscle gap-junction coupling. This theoretical result (Fig. 10A,B) 
agrees with the recent experimental results on unc-9(fcl6) mutants (Liu 
et al, 2006). These mutants have body wall muscle gap junction con- 
ductance severely reduced (but neurotransmission at neuromuscular 
junction unaffected) and velocity much lower than wild type (Liu et al, 
2006) . As velocity and frequency are hnearly related (see Karbowski et 
al, 2006), this implies that in unc-9(fcl6) worms the frequency is also 
strongly reduced, as our theory predicts (Fig. 10A,B). 

• Frequency decreases monotonically with increasing model time constants. 
This theoretical result (Fig. 12) is consistent with the data on unc- 
54 mutants. These mutants have significantly increased time constant 
characterizing the duration of the cross- bridge cycle formation between 
myosin and actin (Moerman and Fire, 1997) and their undulatory fre- 
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quency is much smaller than wild type (Fig. 3A,B; Suppl. Table 
Tl). In our model, Eqs. (1-13), this slowing down of the contraction- 
relaxation cycle corresponds to an increase of the muscle membrane 
time constant t„,. 
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Discussion 



In this paper, we use a combined approach involving genetic perturbations, quanti- 
tative behavioral tracking, and mathematical modeling to uncover the mechanisms of 
C. elegans locomotion. Our circuit model is detailed enough to be biologically realistic 
and yet not excessively complex to obscure the basic mechanisms. We hope this model 
provides a link between molecular/cellular and systems level of understanding. 

Discussion of the main results 

We propose the existence of a mechanosensory feedback in C. elegans that strongly 
affects locomotion. Feedback coupling is implied by our locomotory data on non- 
neuronal mutants (Fig. 4; Suppl. Table Tl), in which neural function is not affected 
directly and yet they exhibit abnormal locomotion as compared to wild type. In par- 
ticular, these non-neuronal mutants have significantly modified frequency and flex of 
undulations, which are indirectly linked to the activation levels of the motor neurons. 
The conclusion about the existence of the mechanosensory feedback is also compatible 
with previous experimental data suggesting that altered mechano-sensory channels can 
profoundly affect locomotory pattern. For example, C. elegans unc-8 mutants, which 
have defective putative mechanosensory channels, display abnormal locomotion (Tav- 
ernarakis et al, 1997). Also recent results on the mechanosensory TRPN channel TRP-4 
show that it is involved in nematode bending control (Li et al, 2006). Experimental data 
indicate that some neurons serve as mechanosensors with specialized sensing channels, 
mostly in the head and tail (Goodman and Schwarz, 2003; Suzuki et al, 2003; Wicks 
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et al, 1996). In the theoretical model of the locomotory circuit, we implemented the 
mechanosensory feedback as a stretch receptor coupling in the head and in the body 
wall. Its role, however, is different in both of these circuits. In the head, the stretch 
feedback participates in the control or generation of oscillations, while in the body wall, 
in addition, it generates the neuromuscular wave by breaking a spatial longitudinal 
symmetry of the system through unidirectional coupling between body wall compo- 
nents. Possible neural candidates for the stretch receptors are: a pair of SAA neurons 
in the head, and VB and DB motor neurons along the body wall, because of their long 
undifferentiated processes (White et al, 1986). 

The remarkable experimental finding (Fig. 3 and Suppl. Table Tl) is that the op- 
posite motions - forward and backward - exhibit spatial gradients of flex and in the case 
of neuronal mutants also weak frequency gradients directed in the same way, i.e., from 
the head to the tail. This finding might suggest that primary oscillations are generated 
in the same location for both types of motion. Where could this location be? Because 
in most cases signals propagating through a biological tissue get attenuated as a result 
of energy dissipation caused by tissue viscosity, the most likely candidate for a primary 
pacemaker (CPG) in C. elegans is the head. Below, we argue that other potential can- 
didates for the primary CPG, i.e. either tail or body wall or some combination, are less 
likely. 

If tail served as a primary oscillator, then according to Fig. 3, this would require 
amplification of the signal as it travels from the tail to the head, which in turn, would 
require large amounts of mechanical energy generated in the body wall to compensate 
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for energy loss. Small number and low linear density of the body wall motor neurons 
(37 forward motion motor neurons (VB, DB, VD, DD) cover ~ 80 % of the body length 
as compared to ~ 200 neurons in the head; see White et al, 1986), and the decaying 
trend of synaptic density toward the tail (Fig. 7B in Chen et al, 2006) are probably 
not sufficient to excessively activate body wall muscles and achieve signal amplification. 
Another possibility is that both head and tail are two separate primary pacemakers with 
two distinct spontaneous frequencies. However, because head and tail have numerous 
recurrent synaptic connections, these two potential oscillators would likely adjust their 
frequencies to a single one, as happens in a system of coupled oscillators (e.g. see 
Strogatz, 1994). Some neuronal mutants show that head and tail oscillate differently, 
which is a strong evidence against such a scenario. It is important to stress that our 
hypothesis about the primary oscillator in the head does not exclude the tail from being 
a secondary oscillator. This potential oscillator could be activated only by an outside 
source, either by head neurons or by a mechanical stimulation (e.g. tail tapping). In 
fact, since the tail has the second largest density of neurons in C. elegans (Chen et al, 
2006), it is potentially capable of a slight signal amphfication, as is evident from flex 
data for some mutants near their tails (signal amplification is possible because a typical 
muscle cell in the tail should get more excitation than a typical muscle cell in the body 
wall due to a higher neural density in the tail) . Finally, there also exists a theoretical 
possibility that the body wall motor neurons serve as a primary CPG. This hypothesis 
is, however, not compatible with our experimental results on posterior gradients of 
body undulations, i.e., one would expect more uniform spatial distributions. Other 
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arguments against this possibility are: the non-zero flex of GABA null mutants (Table 
Tl in the Suppl. Information) and older data on interneuron ablation (Chalfie et al, 
1985). GABA null mutants have functionally disabled ventral-dorsal cross-inhibition at 
the neuromuscular junction along the body wall (such cross-inhibition is necessary for 
body undulations). Despite this deficiency, they have non-zero fiexes, which is a clear 
indication that anti-phase undulations must be imposed from another part of the body. 
By the same token, ablation of the two forward motion interneurons AVB and PVC 
have been shown to virtually abolish (or strongly reduce) body wall undulations and the 
forward motion but still keep intact the head undulations (Chalfie et al, 1985). This 
experiment suggests that forward motion body wall oscillations are not sustainable 
without the presence of connections to the head (which again is in support of our 
hypothesis that the head serves as a primary CPG). We want to stress, however, that 
the body wall motor neurons might serve as a secondary short-term pacemaker when 
stimulated from outside, similar to the tail. 

The most interesting consequences of our theoretical model are the biphasic nature 
of the dependence of frequency on the global strength of synaptic coupling (Fig. 9) and 
on muscle calcium signaling (Fig. 10C,D), as well as the constancy of neuromuscular 
wavelength (Fig. 8C,D). The biphasic relationships, being a direct consequence of the 
non-linearities present in the system, are in agreement with the experimental data on 
GABA signaling (Fig. 5), and acetylcholine signaling (Fig. 6) if combined with a pre- 
viously published kinematic results on hyperactive acetylcholine mutants (Karbowski 
et al. 2006). These results may suggest that wild-type worms have synaptic couphngs 



36 



close to their optimal values. The stability of the neuromuscular wave is determined 
by the ability of the locomotory circuit to maintain the precise phase relationships 
between body components. Our model shows that unidirectional stretch receptor cou- 
pling involving long dendrites of the body wall motor neurons (Fig. 2) provides a 
mechanism for such a coherency, since the phase differences between neighboring mus- 
cles and consequently the wavelength are stable and practically frequency and stretch 
receptor coupling independent (Fig. 8C). The invariance of the wavelength agrees well 
with the experimental data in Table 2 and in Karbowski et al (2006). 

Within our theoretical model, we consider two extreme forms of CPG in the head 
circuit. In model A, the primary oscillations are generated on the level of head in- 
terneurons and then only modulated by stretch receptor feedback, and ultimately by 
the whole locomotory network, which makes oscillations more robust. In model B, os- 
cillations in the head emerge as a result of collective loop interactions between head 
neurons, muscles and stretch receptors, without prior oscillatory activity of neurons. In 
this scenario, the stretch receptor feedback is a critial element in creating oscillations, 
unlike in model A. These two mechanisms can be realized in the same circuit with 
the same topology, only with different parameters characterizing head interneurons in- 
teractions, and worms can potentially use either of these mechanisms by modulating 
their respective parameters (change in only few of them is needed to transit between 
models A and B). The presence of oscillations can be controlled by sensory input that 
can change the level of excitation coming to head interneurons, which is controlled in 
our model by the parameter Ci. In general, both excessive excitation (large positive 
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ci) and excessive inhibition (large negative ci) abolish the oscillations and consequently 
the motion. From a mathematical perspective, oscillations in the circuit are the result 
of non-linear feedback loops, which cause temporal delays in the system and make the 
steady state unstable. This mechanism is quite general and robust (insensitive on the 
precise values of parameters) and appears in many biological and non-biological (e.g. 
engineering) systems (Strogatz, 1994). It is important to stress that both of the above 
potential forms of CPG are indistinguishable on a phenotypic level, if stretch receptor 
coupling is sufficiently strong and slow (Figs. 9-12). In this regime, both models yield 
qualitatively similar predictions regarding dynamic properties of the locomotory circuit 
under different perturbations, and both are consistent with the experimental data. This 
fact suggests that global locomotory output is to a large extent a result of the collec- 
tive dynamics of the whole circuit, not just its local elements, and it does not depend 
on a specific mechanism of oscillation generation in most of the parameter space. In 
our model, the head interneuron circuitry is treated in a simplified manner given the 
complexity of the actual circuitry. This simplified approach, however, captures an im- 
portant feature of interneuron connectivity pattern, namely small-scale feedback loops. 
The model head interneurons X, Y, Z form a minimal circuit (loop) in C. elegans that 
can generate robust oscillations. In general, our results suggest that more experimental 
attention should be paid to the head circuit in the context of locomotion. 

Discussion of additional aspects of this study 

Our theoretical study shows the existence of several "windows of parameters" (usu- 
ally sufficiently wide) for the presence of undamped oscillations along worm's body, 
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which is required for a sustained undulatory motion. Beyond those critical ranges of 
parameters oscillations decay just after a few cycles, with the most pronounced decay 
for posterior parts of the circuit. For example, the global levels of excitation and inhi- 
bition in the network can be neither too small nor too large for the stable oscillations 
(Fig. 7). The weak frequency gradients we observe experimentally for neuronal mu- 
tants with compromised synaptic transmission (Figs. 3, 5 and 6; Suppl. Tables Tl, T2) 
may be directly related to the spatial decay of the waveform along the body, which we 
notice in simulations (Fig. 7A,B lower panels). In particular, our simulations indicate 
that for mutants with weak synaptic coupling the more posterior a given muscle is the 
more damped are its oscillations. It seems that our automated experimental setup has 
problems with detecting low amplitude oscillations (low signal-to-noise ratio) in more 
posterior parts, which manifests itself in weak frequency gradients. 

Similarly, muscle calcium signaling cannot be too large to maintain undulatory lo- 
comotion. Experiments by Maryon et al. (1996, 1998), in which wild-type worms were 
exposed to ryanodine, are consistent with the last prediction. Worms treated with ryan- 
odine, which produces an excessive amount of calcium in muscle cells, are either very 
sluggish or completely paralyzed (Maryon et al, 1996). 

In our model, the activity of AVE and PVC has two other important functions. 
First, to make oscillations in the body wall circuit (imposed from the head via muscles 
and stretch receptors) more robust by releasing excitatory motor neurons (VB and 
DB) from inhibitory influence of DVA (when this inhibition is too strong or a signal 
from AVB and PVC too weak oscillations are not sustained). Second, to start the 
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posterior neuromuscular wave by signaling to VB and DB motor neurons with the help 
of their long posterior dendrites. This wave is a second critical factor for the forward 
locomotion, as is evident from a biomechanical model of worm's undulations (Karbowski 
et al, 2006). Following the Byerly and Russell hypothesis we assume that long dendrites 
contain stretch receptors that provide posterior coupling between spatially separated 
parts of the body wall. Without this couphng, the wave is abolished in both models A 
and B. The involvment of stretch receptors in generation of the wave was also suggested 
before on purely biomechanical grounds (Niebur and Erdos, 1993). It is important to 
stress that forward and backward motion require neuromuscular waves traveling in 
opposite directions, which are correlated with the corresponding opposite directions of 
the motor neuron's dendrites in forward and backward circuits. However, based on the 
above arguments it seems that both types of motion have CPG's located probably in 
the head. 

In the model, spontaneous collective oscillations can be also generated in the body 
wall circuit without the involvement of the head circuit and/or AVB and PVC interneu- 
rons, because of the numerous feedback loops present in the body wall. This mechanism 
was theoretically considered by Bryden and Cohen (2004). We did not, however, pursue 
this possibility for the reasons outlined above. Therefore, the parameters in our model 
are chosen such that the body wall neural circuit does not oscillate in isolation from 
the head circuit, i.e., when AVB and PVC neurons are removed and when there is no 
coupling between muscles in the head and the body wall. In particular, we assumed 
that the excitatory (VB, DB neurons) and inhibitory (VD, DD neurons) motor neurons 
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are intrinsically non-oscillatory. 

Our theoretical results show, for both models A and B, that frequency increases 
weakly with increasing the stretch receptor coupling if it is strong enough (Fig. 11). 
This theoretical result is consistent with data for trp-4 mutants (Li et al, 2006) and 
cuticle sqt-1 mutants. These animals show the opposite behavior, i.e., trp-4 ^-^e faster 
and have larger flex than wild- type (Li et al, 2006), while sqt-1 are slower and have 
smaller flex (Suppl. Table Tl). Decreased (increased) flex angle implies smaller (larger) 
stretching, which in turn suggests weaker (stronger) stretch receptor feedback coupling, 
and according to our theory, lower (higher) frequency. The case of BE109 mutants is 
more complex, since they exhibit higher flex angle and lower frequencies than the wild- 
type worms (Suppl. Table Tl). The likely explanation in this case is that despite higher 
flexes the actual stretch receptor coupling might be smaller than in wild-type. This is 
possible because these mutants are noticeably shorter than wild type (about 20%), 
suggesting a profound structural changes in the cutico-skeletal system that might affect 
the quality of a mechano-sensory signaling. 

C. elegans locomotion seems similar to the locomotion of other systems performing 
undulations. However, there are also some important differences, in particular, other 
undulatory systems such as leech, lamprey, and crayfish have anatomically well defined 
body segments, which serve as local CPG's (Skinner and MuUoney, 1998; Priesen and 
Gang, 2001; Marder et al, 2005). In contrast, C. elegans lacks body segments (White el 
al, 1986) and its body wall seems to be rather passive in comparison to the head, which 
likely serves as CPG. 
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We note that our circuit model allows a direct investigation of the effects of single- 
or multi-neuron ablations (e.g. motor neurons) on locomotion characteristics. By com- 
paring kinematic parameters of C. elegans with ablations to circuit model predictions 
one can in principle gain information on the strength of some synaptic connections and 
refine the model. 

Some of our extrachromosomal transgene constructs involving slo-1 mutants might 
not precisely restore the wild- type functions. However, because our analysis of the 
experimental data is to a large extent qualitative, the possible under- or over-shooting 
to the wild-type level are not critical to our conclusions. 

Conclusions 

Our major conclusions are: (i) based on the anterior-to-posterior gradients of the 
flex for almost all strains, and partly based on the corresponding weak gradients of 
the frequency for neuronal mutants, the most hkely location of CPG for both forward 
and backward locomotion is the head, (ii) stretch receptor feedback associated with 
long unidirectional dendrites in the body wall plays a critical role in neuromuscular 
wave generation and conservation of the wavelength, (iii) there exist optimal values of 
synaptic couplings (both inhibitory and excitatory) for which frequency of undulations 
is maximal, and wild-type worms might be close to that optimal regime, (iv) the precise 
mechanism of oscillation generation within CPG, i.e., whether neural activity or reflex 
activity dominate, is irrelevant and indistinguishable on a phenotypic level, except for 
weak and fast stretch receptor coupling, (v) frequency increases monotonically with 
muscle gap-junction strength, and (vi) frequency decreases monotonically with increas- 
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ing time constants of the circuit elements. Our circuit model can also make predictions: 
(i) frequency should depend in a non-monotonic manner on muscle calcium signaling, 
and (ii) long dendrites of the excitatory motor neurons in the body wall are critical 
for generating neuromuscular wave and thus perturbing them should perturb the wave- 
length. Experiments addressing these problems would either provide additional support 
for our model or cast doubt on some assumptions, thus helping in refining the model. 
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Table 1: Parameters used in simulations of Eqs. (1)-(13). 
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Inhibitory w^z is assumed to be non-GABA-ergic. Labels (A) and (B) refer 
to models A and B. 
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Table 2: Body length normalized wavelength during forward and backward 
locomotion for adult wild-type (WT) C. elegans and mutants. 



Genotype 


Normalized wavelength 






Forward motion 


Backward motion 


a elegans WT (N=62) 


0.61±0.3 


0.64±0.07 


Neuronal Mutants: 






unc-25(el56) (N=4) 


0.68±0.01 


0.69±0.04 


unc-46(el77) (N=4) 


0.71±0.01 


0.71±0.02 


unc-18(e81); cho-l::unc-18::yfp (N=14) 


0.70±0.03 


0.71±0.07 


slo-l(jsll8) (N=18) 


0.62±0.02 


0.64±0.02 


Non-neuronal Mutants: 






egl- 19(n582 ); unc- 1 1 9::egl-19::yfp 






Line A (N=16) 


0.64±0.03 


0.67±0.13 


Line B (N=16) 


0.63±0.03 


0.66±0.07 


unc-68(rll58) (N=16) 


0.62±0.03 


0.62±0.06 
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Figure Legends 
Fig. 1 

The general approach taken in this paper. (A) Wc perturb (manipulate) genes whose 
products affect some properties of different classes of neurons and other non-neuronal 
parameters. These altered properties, in turn, affect locomotory characteristics. The 
locomotory pattern may affect activities of neurons via mechanosensory feedback. (B) 
Actual image of a worm with superimposed reference points measuring worm's coordi- 
nates. Based on a spatio-temporal pattern of each of these points we calculated motion 
characteristics. The maximal angle over time between two solid lines (green and red 
on-line) determines the half of a local flex angle, which quantifies a deviation from a 
straight line of two neighboring reference body sections during motion. Dashed arrow 
(blue on-line) indicates the direction of motion). Adapted from Cronin et al. (2005). 

Fig. 2 

Schematic diagram of the circuit for forward locomotion in the nematode C. elegans. 
Excitatory synapses are denoted by arrows (— >), inhibitory synapses by bars (H), and 

electric gap junctions by dashed hues ( ). (A) The large scale view of the circuit. 

Head neurons excite and inhibit head muscles (ventral M^^.^ and dorsal Mm), and ad- 
ditionally receive feedback from head stretch receptors (ventral Shv and dorsal Shd)- 
The signal from the head neurons is transmitted to the body wall circuit through two 
interneurons AVB and PVC (with very long dendrites spanning the whole body length), 
and through the gap junction coupling between head and body wall (My and Ma) mus- 
cles. Body wall excitatory motor neurons (ventral VB and dorsal DB) use acetylchohne 
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for synaptic signaling. These neurons have elongated processes directed posteriorly, 
which, we assume, contain stretch receptors (S). Stretch receptors detect relative dif- 
ferences in the activities of the corresponding ventral and dorsal muscles. Body wall 
inhibitory motor neurons (ventral VD and dorsal DD) use GABA for synaptic signaling. 
These neurons act as cross-inhibitors, since their only input is from the opposite side 
excitatory motor neurons. Excitatory motor neurons as well as inhibitory motor neu- 
rons are connected within each group by gap junctions (not shown). Also, excitatory 
motor neurons (VB and DB) receive synaptic input from DVA interneuron, which we 
assume is inhibitory. (B) Schematic diagram of the head CPG. The head interneurons 
X (corresponding to SAA neuron) and Y are excitatory and, we assume, signal with 
acetylcholine, while Z is inhibitory and, we assume, signals without GABA (GABA is 
expressed only in one head interneuron RIS and its ablation does not influence pheno- 
type; see Schuske et al, 2004). Ventral and dorsal groups of these interneurons effectively 
inhibit each other such that oscillatory rhythm in the head alternates between ventral 
and dorsal sides. Head excitatory motor neurons {Ehy and Ehd corresponding, e.g., to 
pairs of SMB and SMD neurons) use acetylcholine, and head inhibitory motor neurons 
{Ihv and Ihd corresponding to a pair of RME neurons) use GABA for synaptic signaling. 
As in the body wall, stretch receptors in the head {Shv and Sm) detect relative differ- 
ences in the ventral and dorsal activities of head muscles {M^v and M^d}- The primary 
oscillations can be generated either within the X, F, Z loop or with a stretch recep- 
tor feedback within the X, Y, E, M, S loop. Changes in the level of global excitation 
(acetylcholine) and inhibition (GABA) affect frequency of these oscillations. 
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Fig. 3 

Experimental examples of distributions of the undulatory frequency and bending flex 
angle along the worm's body for forward and backward motion. (A), (B) Frequency is 
relatively uniform during forward and backward motion for wild type (solid blue line) 
and non-neuronal mutants unc-68(rll58) (dashed-dotted green line) and unc-54(s74) 
(dotted black hue), but it decays weakly towards the tail for GABA defective unc- 
25(el56) mutant in forward locomotion (dashed red hne). (C), (D) Flex angle decays 
markedly posteriorly for all these animals, in some cases, with an increase towards the 
tail (notably for unc-54(s74 ))■ The correlation coefficient r between the flex and position 
from the head has the following values. For forward motion: r = —0.97 (p = 0.000) 
for wild type; r = -0.86 {p = 0.001) for unc-25(el56); r = -0.83 {p = 0.002) for 
unc-68(rll58); r — —0.71 (p = 0.014) for unc-54(s74)- Similarly, for the backward 
motion: r = -0.89 (p = 0.000) for wild type; r = -0.62 (p = 0.042) for unc-25(el56)] 
r = -0.72 (p = 0.013) for unc-68(rll58); r = 0.41 (p = 0.208) for unc-54(s74). In 
all cases but one (for unc-54 the backward motion r is positive but not significant) 
correlation are strong and negative, which implies that there is a significant decreasing 
trend of flex with position along the body. The line convention is the same as in (A), 
(B). Here and in other figures the worm's head position corresponds to the point 1 on 
the horizontal axis and its tail position to the point 11. The error bars in all figures are 
population standard deviations. 
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Fig. 4 

Locomotion characteristics of the mutants with affected non-neuronal functions and 
wild- type C. elegans (forward motion). (A) Spatial distribution of the frequency of un- 
dulations for cuticle mutants sqt-l(scl01 ) (dashed hne (red on-hne), = 13) and BE109 
(dashed-dotted line (green on-line), A" = 9), and wild- type (solid line (blue on-line), 
A" = 62) worms. The error bars are population standard deviations. (B) Snapshots of 
worm postures: (a) wild type, (b) BE109, (c) sqt-l(sclOl), (d) unc-54(s95). All of the 
mutants move on average slower than wild type, but their body posture and locomotory 
pattern can be either similar to wild type (panel (c)) or significantly different (panels 
(b) and (d)). 

Fig. 5 

Experimental non-monotonic (biphasic) dependence of the undulatory frequency on 
synaptic GABA signaling (forward motion). (A) Loss-of-function GABA mutants unc- 
25(el56) (dashed red line, A^ = 4), unc-46(el77) (dashed-dotted green hne, A^ = 4), 
and unc-18(e81); cho-l::unc-18::yfp (dotted black line. A" = 14) move with significantly 
lower frequencies than wild type (solid blue line) (non-parametric sign test for paired 
average frequencies along the body, p < 0.05). (B) Mutants slo-l(jsll8); Pacr-2::slo-l 
(dashed-dotted black lines representing 4 independent lines A, B, C, D) with increased 
GABA signaling and presumably wild-type levels of acetylcholine signaling move with 
significantly lower frequencies than wild type (solid blue hne) (sign test, p < 0.05). 
(C) Similar to panel B: mutants slo-l(jsll8); Pacr-2::slo-l (dashed-dotted black hne) 
are also slower than two types of animals with presumably recovered wild-type neu- 



54 



rotransmission in both GABA and cholinergic neurons: slo-l(jsll8); Pacr-2::slo-l; 
Punc-25::slo-l (sohd red hne, N = 21) and slo-l(jsll8); Psnh-l::slo-l (sohd green 
hne, N = 20) (sign test, p < 0.05). 

Fig. 6 

Experimental dependence of the undulatory frequency on synaptic acetylcholine signal- 
ing (forward motion). (A) In a background with elevated GABA signahng, slo-l(jsll8); 
Pacr-2::slo-l worms with presumably wild-type levels of acetylcholine move with either 
slightly higher (lines A {N — 19) and B {N — 18); dashed red line) or shghtly lower 
(lines C {N = 18) and D (A^ = 17); dashed green line) frequencies than the control 
slo-l(jsll8); Pacr-2::gfp (solid blue line, N = 19), which has both types of neurotrans- 
mission increased above the wild-type level. (B) Mutants with increased acetylcholine 
and presumably rescued GABA signaling towards wild-type level, i.e., slo-l(jsll8); 
Punc-25::slo-l (dashed red line; N = IS (line A), N = IQ (line B)), slo-l(jsll8); Punc- 
17::slo-l (dashed-dotted green line; AT = 16 (hne A), A^ = 18 (hne B), A^ = 18 (hne C)), 
and slo-l(]sll8); Pcho-l::slo-l (dotted black line; N = 1% (line A), AT = 18 (line B), 
A" = 17 (line C)), have slightly lower frequencies that the wild-type worms (solid blue 
line) (nonparametric sign-test for paired average frequencies along the body, p < 0.05). 

Fig. 7 

Oscillations generated in the circuit. (A) Model A. Upper panel: ventral muscle ac- 
tivities along the circuit (head Mhv - solid line (blue on-line), M„3 - dashed line (red 
on-line), M^q - dashed-dotted line (green on-line)). Middle panel: corresponding dorsal 
muscle activities {M^^, M^, M^e) along the circuit. Note alternate corresponding os- 
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cillations in both these panels {qex — 4.0, qin — 3.0). Lower panel: damped oscillations 
along the circuit (the same line convention as in the upper panel). For more posterior 
circuit elements oscillations become unstable only after a few cycles when either global 
excitation and inhibition are too small or muscle gap-junction coupling is too small 
{qex — 2.0, qin = 1.0). This damping effect may explain the experimentally observed 
weak spatial decay of frequency for neuronal mutants with weak synaptic couplings 
(Suppl. Tables Tl, T2). (B) The same for model B. (C) Oscillations in the circuit ele- 
ments of model A. Upper panel: activities of head ventral X (solid line (blue on-line)) 
and Z (dashed line (red on-line)) neurons. Middle panel: activities of two interneurons 
AVB (solid line (blue on-line)) and PVC (dashed hne (red on-line)). Note weak high 
frequency oscillations in both interneurons in contrast to other elements in the circuit. 
Lower panel: activities of ventral body wall motor neurons (no 2) excitatory (solid line 
(blue on-line)) and inhibitory (dashed line (red on-line)), and stretch receptors (dashed- 
dotted line (green on-line)) {qex = 3.0, gj„ = 2.0). Note that almost all elements in the 
circuit oscillate with the same collective frequency, except the two interneurons, which 
exhibit weak high frequency oscillations (for some model parameters these oscillations 
are practically invisible). (D) The same for model B. 

Fig. 8 

Neuromuscular wave generated in the circuit. The wave propagates either towards 
the tail (A) for posteriorly directed stretch receptor coupling, or towards the head 
(B) for anteriorly directed coupling. The former corresponds to the forward motion 
circuit, whereas the latter to the backward motion circuit. Wave does not propagate 
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when a unidirectional stretch receptor feedback couphng is absent. (C) The phase 
lag of neighboring muscle activities is to a large extent both position and frequency 
independent. Also the phase lag between the head and the tail is essentially frequency 
independent: 0.97 x (27r) for 0.70 Hz, 1.03 x {2tt) for 0.57 Hz, and 0.91 x {2tt) for 0.46 Hz, 
which suggests frequency independence of the wavelength. Solid line (blue on-line) with 
diamonds corresponds to the frequency 0.70 Hz, dashed line (red on-line) with circles to 
0.57 Hz, and dashed-dotted hne (green on-line) with crosses to 0.46 Hz. (D) The phase 
lag of intermuscle activities (between no. 3 and 2) decreases with increasing the stretch 
receptor feedback coupling, although weakly if the coupling is sufficiently strong for 
both models A (sohd line (blue on-line)) and B (dashed hne (red on-hne)). Parameters 
used in (A) and (B): = 4.0, gj„ = 2.0, and gm = 0.4 (model A). Parameters used in 
(C) (model B only) and (D) (both models): Qex — 3.0, Qin — 2.0. 

Fig. 9 

Theoretical dependence of the frequency of oscillations in the circuit on the global 
synaptic inhibition gj„ (A, B) and global synaptic excitation q^^ (C, D). Both models A 
and B give a qualitatively similar biphasic dependence of frequency on inhibition and 
excitation. (A), (B) For very weak inhibition oscillations are unstable. (C), (D) For 
both too weak and too strong excitation oscillations are unstable. (A), (B) The solid 
line (blue on-line) corresponds to Qex = 4.0 and dashed line (red on-line) to Qex = 2.0. 
(C), (D) The sohd hne (blue on-hne) corresponds to Qin — 5.0 and dashed line (red 
on-line) to q^n = 1.0. These biphasic relationships can be explained by the changes 
in temporal characteristics of muscle voltage in response to changes in the levels of 
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inhibition (E) and excitation (F). Parameters used: (E) — 4.0 and qin — 0.5 (upper 
panel), Qin = 4.5 (middle panel), Qin = 35.0 (lower panel) in model B; (F) q^n = 5.0 and 
qex = 1.2 (upper panel), q^n — 3.8 (middle panel), q^n — 6.6 (lower panel) in model A. 

Fig. 10 

Theoretical dependence of the frequency of oscillations in the circuit on the muscle 
characteristics. Models A and B give qualitatively similar behaviors of frequency on 
muscle gap junction coupling (A, B) and muscle calcium signaling strength Wjnm (C, 
D). (A), (B) Monotonic increase of frequency with increasing gap junction coupling. 
Solid hne (blue on-line) corresponds to qex = 2.0, q^ = 2.0 and dashed line (red on- 
line) to qex = 2.0, qin — 4.0. (C), (D) Biphasic or triphasic dependence on frequency 
on calcium signaling. The curve's shape depends on a relative strengths of the global 
excitation and inhibition. For very large calcium signaling levels oscillations in the 
circuit become unstable. Solid line (blue on-line) corresponds to qex = 4.0, q^n = 4.0 
and dashed line (red on-line) to qex — 1-5, qin — 7.0. 

Fig. 11 

Theoretical dependence of the frequency of oscillations in the circuit on the stretch 
receptor coupling parameter Wes- The behavior of model A (A) differs from model B 
(B) only in the weak coupling regime. In this regime oscillations in model B are unstable, 
while model A shows unmodulated high frequency oscillations of head interneurons. In 
both figures, the solid line (blue on-line) corresponds to qex = 2.0, qin = 2.0 and the 
dashed line (red on-line) to qex — 4.0, qin — 1.0. 
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Fig. 12 

Theoretical dependence of the frequency of oscillations on the muscle and stretch re- 
ceptor time constants. (A) Decrease of the oscillatory frequency with increasing the 
muscle membrane time constant in model A (solid line (blue on-line)) and model B 
(dashed line (red on-line)). (B) Decrease of the oscillatory frequency with increasing 
the stretch receptor time constant in model A (solid line (blue on-line)) and model 
B (dashed hne (red on-hne)). Note that the frequency is much larger in model B for 
small values of Tm and Ts- In both figures, we used q^x = 3.0 and qin = 2.0 
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